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Abstract. We present the perturbative parts of the structure functions F2 and F£ for 
a gluon target having nonzero transverse momentum squared at order as- The results 
of the double convolution (with respect to the Bjorken variable x and the transverse 
momentum) of the perturbative part and the unintegrated gluon densities are compared 
with HERA experimental data for F2 and Fl at low x values and with predictions of 
other approaches. The contribution from Fl structure function ranges 10-^30% of that 
of F2 at the HERA kinematical range. 

1 Introduction 

The basic information on the internal structure of nucleons is extracted from 
the process of deep inelastic (lepton-hadron) scattering (DIS). Its differential 
cross-section has the form: 



where F2{x,Q'^) and Fl{x,Q'^) are the transverse and longitudinal structure 
functions (SF) , respectively, and are the photon and the hadron 4-momentums 
and X = Q'^/(2pq) with = _q2 > g. 

In the lecture we will study only the charm part F| of the transverse SF F2 
and the longitudinal SF F^. 

The study is related with the fact that recently there have been important 
new data on the charm SF F2, of the proton from the HI [1] and ZEUS [2] 
Collaborations at HERA, which have probed the small- a; region down to .t = 8 x 
10~^ and X — 2x 10^'^, respectively. At these values of x, the charm contribution 
to the total proton SF, F2, is found to be around 25%, which is a considerably 
larger fraction than that found by the European Muon Collaboration (EMC) 
at CERN [3] at larger x, where it was only ~ 1% of F2. Extensive theoretical 
analyses in recent years have generally served to confirm that the F2 data can 
be described through perturbative generation of charm within QCD (see, for 
example, the review in Ref. [4] and references therein). 

The second object of our study is the longitudinal SF F]^{x, Q^)- It is a very 
sensitive QCD characteristic because it is equal to zero in the parton model 



dxdy 



xQ^ 
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with spin— 1/2 partons. Unfortunately, essentially at small values of x, the ex- 
perimental extraction of Fl data requires a rather cumbersome procedure (see 
[5,6], for example). Moreover, perturbative QCD leads to some controversial re- 
sults in the case of SF Fl. The next- to-leading order (NLO) corrections to the 
longitudinal coefficient function, which are large and negative at small x [7,8], 
need a resummation procedure ^ that leads to a coupling constant scale higher 
essentially than (see [10,8,11]) ^. 

Recently there have been important new data [17]-[22] of the longitudinal 
SF Fl, which have probed the small- a; region down to a; ~ 10"-^. Moreover, the 
SF Fl can be related at small x with SF F2 and the derivative dF2/dlnQ^ (see 
[23]). In this way most precise predictions based on data of F2 and dF2/d\nQ'^ 
(see [19] and references therein) can be obtained for Fl- These predictions can 
be considered as indirect 'experimental data' for Fl- 

We note, that perhaps more relevant analyses of the HERA data, where 
the x values are quite small, are those based on BFKL dynamics [16], because 
the leading ln(l/a:;) contributions are summed. The basic dynamical quantity 
in the BFKL approach is the unintegrated gluon distribution (pg{x,k'^) ^ (fg 
is the (integrated) gluon distribution multiplied by x and k± is the transverse 
momentum) 

fgix.Q'') = fg{x,Ql) + / dkl^g{x,kl), (1) 

Jqi 

which satisfies the BFKL equation. The integral is divergent at the lower limit 
and it leads to the necessity to use the difference fg{x, Q^) — fg{x, Qq) with some 
nonzero Qq. 

Then, in the BFKL-like approach (hereafter the fct-factorization approach 
[24,25] is used) the SF i^|^(a;, Q^) and Fl{x, Q"^) are driven at small x by gluons 
and are related in the following way to the unintegrated distribution ipg{x,kj_): 

F^lM^) = f'^ jdklel Cl^{z,Q\mlkl) <Pg{x/z,kl), (2) 
Fl{x,Q^) -Fl{x,Q^) 

dkl J2 elCl^{z,Q\0,kl)^g{x/z,kl), (3) 

i=u^d,s 

where e, is the charge of the i-flavor quark. 

^ Without a resummation the NLO approximation to Fl can be negative at low x and 

quite low values (see [8,9]). 
^ Note that at low x a similar property has been observed also in the approaches [12- 

14] (see recent review [15] and discussions therein), which based on Balitsky-Fadin- 

Kuraev-Lipatov (BFKL) dynamics [16], where the leading ln(l/a;) contributions are 

summed. 

^ Hereafter is the gluon 4-momentum. 
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Fig. 1. The diagrams contributing to T^j^ for a gluon target (the upper and down wave 
lines on the diagram (a) correspond to photons and gluons, respectively). They should 
be multiplied by a factor of 2 because of the opposite direction of the fermion loop. 
The diagram (a) should be also doubled because of crossing symmetry. 



The functions l{x,Q'^ ,m1,k\) and C£(a;, Q^, 0, A;^^) may be regarded as 
the structure functions of the off-shell gluons with virtuality k\ (hereafter we call 
them as hard structure functions They arc described by the quark box (and 
crossed box) diagram contribution to the photon-gluon interaction (see Fig.l.) 

The purpose of the paper is to present the results of [26] for these hard 
SF C|^(x,Q^,m^,A:5.) and to analyze experimental data for F|(a;, Q^) and 
Fhi^jQ"^) by applying Eq. (2) with different sets of unintegrated gluon den- 
sities (see Ref. [27,15]) and to give predictions for the longitudinal charm SF 
F£(x,g2). 

It is instructive to note that the results should be similar to those of the 
photon-photon scattering process. The corresponding QED contributions have 
been calculated many years ago in Ref. [28] (see also the beautiful review in 
Ref. [29]). Our results have been calculated independently in [26] (based on 
approaches of [30,31]) and they are in full agreement with [28]. However, we 
hope that our formulas which are given in a simpler form could be useful for 
others. 



* This notation reflects the fact that SF F2_^ and Fl connect with the functions C| ^ 
and C£ at the same form as cross-sections connect with hard ones (see [24,25]). 
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2 Hard structure functions 



The gluon polarization tensor (hereafter the indices a and /? are connected 
with gluons), which gives the main contribution at high energy Hmit, has the 
form: 



pa/9 
^BFKL 



Ua jL/3 



-A;2 



1 1 

2^ 



(4) 



where = 1 - 4bx'^, b = -P/Q'^ = k'i/Q'^ > 0. 

Contracting the corresponding photon projectors, we have: 



f (1) + _ f (2) 

2/32 



(5) 



/(I) = i [^2^^(l) _ 352.2/(1)1 y(2) ^ ^ 



where the as{Q^) = aa{Q'^)/{A-K), a = /Q"^ and 



/(I) = -2/3 



1 - ^1 - 2a;(l + 6 - 2a) [1 - a;(l + 6 + 2a)]^ /i 



+ (2a-6)(l-2a)a;V2 



(6) 



= 8a; /3 



(1 - (1 + b)x) - 2x(^x{l - (1 + b)x){l + b-2a) + a^^^ /i 



+ 6a;2(l-(l + 6)ar)(2a-6)/2 



(7) 



l-a;(l + 6) 



2[x{l-x{l + b)){l + b-2a) + a~p'' fi 



-x{l-x{l + b)){l-2a) f2 



= 4 /3 (1 - (1 + b)xf 



2 - (1 + 26a;2) /i - bx^ /2 



(8) 



(9) 



with 



Aax 



{i-{i + b)xy 



1 1 + /3/3 
/i = — In ^, 



h = 



-4 



1 - /32/32 



For the important regimes when k"^ = 0, ml — Q and/or = 0, the results 
coincide with ones in Refs. [24,32]. Notice that our results in Eq. (5) should also 
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agree with those in Ref. [33] but the direct comparison is quite difficult because 

the structure of their results is quite cunibcrsonic (sec Appendix A in Ref. [33]). 
We have found numerical agreement in the case of F2{x, Q^) for several types of 
unintegrated gluon distributions (see Fig. 4 in [26]). 

3 Relations between Fl, F2 and derivation of F2 in the 
collinear approximation 

Another information about the SF can be obtained in the collinear ap- 
proximation (i.e. when fc^ = 0) in the following way (see also our study [34]). 

In the framework of perturbative QCD, there is the possibility to connect Fl 
to F2 and dF2 /rf In due the fact that at small x the DIS structure functions 
depend only on two independent functions: the gluon distribution and the singlet 
quark one (the nonsinglet quark density is negligible at small x), which in turn 
can be expressed in terms of measurable F2 and dF2/dlnQ'^. 

In this way, by analogy with the case of the gluon distribution function (see 
[35,36] and references therein), the behavior of Fi(a;, Q"^) has been studied in [23], 
using the HERA data [37,38] and the method [39] ^ of replacement of the Mellin 
convolution by ordinary products. Thus, the small x behavior of _FY(x, Q^) can 
be extracted directly from the measured values of F2{x,Q'^) and its derivative 
without a cumbersome procedure (see [5,6]). These extracted values of Fl may 
be well considered as new small x 'experimental data' of Fl- The relations can 
be violated by nonperturbative corrections like higher twists (see [41,42]), which 
can be large in the case of Fl [43,44]. 

Because fcy-factorization approach is one of popular perturbative approaches 
used at small x, it is very useful to compare its predictions with the results of [23] 
based on the relations between SF Fl{x, Q^), F2{x, Q^) and dF^ix, (5^)/dlnQ^. 

The fcy-factorization approach is strongly connected to the Regge-likc behav- 
ior of parton distributions. So, we restrict our investigations to SF and parton 
distributions with the following form (hereafter a = q,g): 

fa{x,Q^)^F2{x,Q^)r.x-'^'^'^ (10) 

Note that really the slopes of the sea quark and gluon distributions: 5q and 
Sg, respectively, and the slope Sf2 of F2 are sligtly different. The slopes have 
a familiar property 6q < 6f2 < 5g (see Refs. [45]-[50] and references therein). 
We will neglect, however, this difference and use in our investigations the ex- 
perimental values of 5{Q^) = Sf2{Q^) extracted by HI Collaboration ^ (see [45] 

® The method is based on previous investigations [5,40]. 

^ Now the preliminary ZEUS data for the slope dlnF2/dln{l/x) arc available as some 
points on Figs. 8 and 9 in Ref. [51]. Moreover, the new preliminary HI points have 
been presented on the Workshop DIS2002 (see [21]). Both the new points axe shown 
quite similar properties to compare with HI data [45]. Unfortunately, tables of the 
ZEUS data and the new HI data are unavailable yet and, so, the points cannot be 
used here. 
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and references therein). We note that the Q^-dependence is in very good agree- 
ment with pcrturbativc QCD at > 2 GcV^ (sec [52]). Moreover, the vahies 
of the slope 6{Q^) are in agreement with recent phenomenological studies (see, 
for example, [12]) incorporating the NLO corrections [53] (see also [54]) in the 
framework of BFKL approach. 

Thus, assuming the Regge-like behavior (10) for the gluon distribution and 
F2{x,Q^) at ^ 1 and using the NLO approximation for coUinear coeffi- 
cient functions and anomalous dimensions of Wilson operators, the following 
approximate results for Fl{x, Q^) has been obtained in [34]: 



r{l + 6)m) 



(l + 30a.(Q2)[i/5_^p,(5)]) 



dF2{xm,Q^) 



+ ^ P2i6)a.iQ^)F2ixaS),Q^) 



rflnQ2 

+ 0{al,asX,x'^), 



where 



r{S) 



AS 



2 + d + 



pi (5) = 1 + 5 + ^2/4, P2{S) 



2.39(5 + 2.69(5^ 



Here 



1 _ 1 
1~S 



_ r{i^5)r{i + u) s 
r{i-s + u) 



(11) 



(12) 



(13) 



where F is the Rimmanian i^-fiuiction. The value of comes (see [55]) from 
asymptotics of parton distributions fa{x) at a; ^ 1: /« ~ (1 — a;)''", and i/ « 4 
from quark account rules [56]. 

Note that the 1/S coincides approximately with 1/6 when 5^0 and a; — > 0. 
However, at 5 — > 0, the value of 1/d is not singular: 



l^\n(^^^ - 'F{l + iy)-^{l) , 
where •Z'-function is the logarithmic derivation of the F-function. 



(14) 



4 Comparison with and Fl experimental data and 
predictions for F^ 

With the help of the results obtained in the previous sections we have ana- 
lyzed HERA data [2] for SF F| from ZEUS Collaboration and the data for SF 
Fl mostly from HI Collaboration. 

In our formula (13) we are mostly interested in gluons, so we take z/ = « 4 below. 
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Notice that in Ref. [24] the fcj^-integral in the r.h.s. of Eq.(2) has been evalu- 
ated using the BFKL results for the Mellin transform of the unintcgratcd gluon 
distribution and the Wilson coefficient functions have been calculated for the 
full perturbative series at asymptotically small x values. Since we would like 
to analyze experimental data for F2 and on a quite broad range of small x 
values, we need a parameterization of unintegrated gluon distribution. 

To study F2 we consider two different parametrizations for the unintegrated 
gluon distribution (see [27]): the Ryskin-Shabelski (RS) one [57] and the Blum- 
lein one [58]. 

We see in Fig. 2 that at large {Q^ > 10 GeV^) the SF F| obtained in 

the kx factorization approach is higher than the SF obtained in pure pertur- 
bative QCD with the GRV [59] gluon density at the LO approximation (see 
curve 1) and has a more rapid growth in comparison with perturbative QCD 
results, especially at ^ 130 GeV^ [60]. For < iq geV^ the predictions 
from perturbative QCD and those based on the kx factorization approach are 
very similar * and show a disagreement with data below Q"^ = 7 GeV^ ^. Un- 
fortunately the available experimental data do not permit yet to distinguish the 
kr factorization effects from those due to boundary conditions [57]. 

The results for the SF Fl obtained in perturbative QCD and from the kr 
factorization approach are quite similar to the F| case discussed above. The 
ratio R"^ = i^£/i^| is shown in Fig. 3. We see that R'^ « 0.1-^0.3 in a wide region 
of Q"^. We would like to note that these values of contradict the estimation 
obtained in Ref. [2]. So, the effect of the large values should be considered in 
the extraction of F2 from the corresponding differential cross-section in future 
more precise measurements. 

For the ratio R'^ we found quite flat .x-behavior at low x in the low 
region (see Fig. 3), where approaches based on perturbative QCD and on kx 
factorization give similar predictions (see Fig 2). It is in agreement with the 
corresponding behaviour of the ratio R = Fl/{F2 — Fl) (see Ref. [23]) at quite 
large values of 6 {6 > 0.2 — 0.3). The low x rise of R'^ at high disagrees 
with early calculations [23] in the framework of perturbative QCD. It could be 
due to the small x resummation, which is important at high (see Fig 2). We 
plan to study this effect in future. 

In Fig. 4 we show the SF Fl as a function x for different values of in 
comparison with HI experimental sets: the old one [17] (black triangles), the last 
year set [19] (black squares) and the new preliminary data [22] (black circles) and 
also with NMC [61] (white triangles), CCFR [62] (white circles) and BCDMS [63] 

® This fact is also due to the quite large value of Qo = 4 GeV^ chosen here. 

^ A similar disagreement with data at < 2 GeV^ has been observed for the complete 
structure function F2 (see, for example, the discussion in Ref. [49] and reference 
therein) . We note that the insertion of higher-twist corrections in the framework of 
usual perturbative QCD improves the agreement with data (see Ref [41]) at quite 
low values of Q^. 

The behaviour is in agreement with previous studies [55,23]. Note that at small 
values of S, i.e. when ~ Const, R has the strong negative NLO corrections (see 
Refs. [7,8]) and tends to zero at x — » after some resummation done in [8]. 
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Fig. 2. The structure function F2{x, Q^) as a function of x for different values of 

compared to ZEUS data [2]. Curves 1, 2 and 3 correspond, respectively, to pure pertur- 
bative QCD with the GRV [59] gluon density at the leading order approximation and 
Ut factorization with RS [57] and Blumlein (at Qo = 4 GeV^) [58] parameterizations 
of unintegrated gluon distribution. 
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Fig. 3. The ratio -R" = Fl{x,Q^)/Fi as a function of x for different values of Q^. 
Curves 1, 2 and 3 are as in Fig. 2. 
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data (white squares). For comparison with these data we present the results of 

the calculation with three different parameterizations for the unintegrated gluon 
distribution ipg{x,k]_,Ql) at Qq = 4 GeV^: Kwiecinski-Martin-Stasto (KMS) 
[64], Blumlein [58] and Golec-Biernat and Wusthoff (GBW) [65]. 

We also add the 'experimental data' obtained using the relation between 
SF Fl{x,Q'^), F2{x,Q'^) and dF2{x,Q^)/d\nQ^ (see Section 3) as black stars. 
Because the corresponding data for F2{x, Q'^) and dF2{x,Q'^)/dlnQ^ are essen- 
tially more precise (sec [19]) to compare with the preliminary data [22] for F^, 
the 'experimental data' have strongly suppressed uncertainties. As it is shown in 
Fig. 4 there is very good agreement between the new preliminary data [22], the 
'experimental data' and predictions from perturbative QCD and fcr-factorization 
approach. 

The differences observed between the curves 2, 3 and 4 are due to the differ- 
ent behavior of the unintegrated gluon distribution as a function of x and k±. 
We see that the SF Fl obtained in the /cT-factorization approach with KMS 
and Blumlein parameterizations is close each other and higher than the SF ob- 
tained in the pure perturbative QCD with the GRV gluon density at the leading 
order approximation. Otherwise, the /cT-factorization approach with GBW pa- 
rameterization is very close to pure QCD predictions: it should be so because 
GBW model has deviations from perturbative QCD only at quite low values. 
Thus, the predictions from perturbative QCD and those based on kr factoriza- 
tion approach are in agreement with each other and with all data within modern 
experimental uncertainties. So, a possible high values of high-twist corrections 
to SF Fl predicted in [44] can be important only at low values: < Qq = i 
GeV^. 

There are several other popular parameterizations (see, for example, Kimber- 
Martin-Ryskin (KMR) [66] and Jung-Salam (JS) [67]), which are not used in our 
study mostly because of technical difficulties. Note that all above parameteri- 
zations give quite similar results excepting, perhaps, the contributions from the 
small /c^-range: k^ < 1 GeV^ (see Ref. [15] and references therein). Because we 
use Qq = 4 GeV^ in the study of SF F2 and Fl, our results depend very slightly 
on the the small fc^-range of the parameterizations. In the case RS, Blumlein, 
GBW and KMS sets this observation is supported below by our results and we 
expect that the application of KMR and JS sets should not strongly change our 
results. 

5 Conclusions 

We have presented the results for the perturbative parts of the SF and 
F£ for a gluon target having nonzero momentum squared, in the process of 
photon-gluon fusion. 

We have applied the results in the framework of kr factorization approach to 
the analysis of present data for the SF F^ and for the charm contribution to F2 
and we have given the predictions for F£ . The analysis has been performed with 
several parameterizations of unintegrated gluon distributions, for comparison. 
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Fig. 4. The structure function Fl{x,Q^) as a function of x for different values of 
compared to experimental data. The HI data: the first 1997 ones [17], the new 
2001 set [19] and the preliminary data [22] are shown as black triangles, squares and 
circles, respectively. The NM [61], CCFR [62] and BCDMS [63] data are shown as white 
triangles, circles and squares, respectively. The 'experimental data' are added as black 
stars. Curves 1, 2, 3 and 4 correspond to SF obtained in the perturbative QCD with 
the GRV [59] gluon density at the leading order approximation and to SF obtained 
in the kr factorization approach with KMS [64], Blumlein (at Qq = 4 GeV^) [58] and 
GBW [65] parametrizations of unintegrated gluon distribution. 
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For SF F2, we have found good agreement of our results with experimental 
HERA data, except at low {Q^ < 7 GeV^). We have obtained also quite large 
contribution of the SF at low x and high {Q^ > 30 GeV^) and this effect 
should be considered for the extraction of from the corresponding differential 
cross-section in future more precise measurements. 

We would like to note the good agreement between our results for and 
the ones obtained in Ref. [68] by Monte-Carlo studies. Moreover, we have also 
good agreement with fits of HI and ZEUS data for F| (see recent reviews in 
Ref. [69] and references therein) based on perturbative QCD calculations. But 
unlike to these fits, our analysis uses universal unintegrated gluon distribution, 
which gives in the simplest way the main contribution to the cross-section in the 
high-energy limit. 

For SF Fl , we have found good agreement between all existing experimental 

data, the predictions for Fl obtained from the relation between SF Fl{x,Q^), 
F2{x,Q-^) and dF2{x,Q^)/dlnQ'^ and the results obtained in the framework of 
perturbative QCD and ones based on fcr-factorization approach with the three 
different parameterizations of unintegrated gluon distributions. 

It could be also very useful to evaluate the complete F2 itself and the deriva- 
tives of F2 with respect to the logarithms of 1/a; and Q"^ with our expressions 
using the unintegrated gluons. 

The consideration of the SF F2 in the framework of the leading-twist approx- 
imation of perturbative QCD (i.e. for "pure" perturbative QCD) leads to very 
good agreement (see Ref. [49] and references therein) with HERA data at low 
X and > 1.5 GeV^. The agreement improves at lower when higher twist 
terms are taken into account [41]. As it has been studied in Refs. [49,41], the SF 
F2 at low is sensitive to the small-x behavior of quark distributions. Thus, 
our future analysis of F2 in a broader range in the framework of kx factoriza- 
tion should require the incorporation of parametrizations for unintegrated quark 
densities, introduced recently (see Ref. [66] and references therein). 
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